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ABSTRACT 

We investigate the use of relativistically blurred photoionized disc models on an XMM- 
Newton observation of the Narrow Line Seyfert 1 galaxy PG 1404+226. The model is 
designed to reproduce the radiation from the inner accretion disc around a Kerr black 
hole, and is more successful at fitting the spectrum than models based on a thermal 
soft excess. The source varies strongly over the course of the observation, and the disc 
model works over all observed flux states. We conclude that it is a useful tool in the 
study of certain quasars. 

Key words: accretion, accretion discs - galaxies: active - galaxies: Seyfert - galaxies, 
individual: PG 1404+226 - X-rays: galaxies. 



1 INTRODUCTION 

Narrow Line Seyfert 1 (NLS1) galaxies are a class of Active 
Galactic Nuclei (AGN) first studied by Osterbrock & Pogge 
(1985). The class is defined by its optical properties: the 
Balmer lines such as H/3 have a FWHM < 2000 km/s, the 
ratio of the [O III] A5007 A forbidden line to the H/3 lines 
is < 3, and emission lines from Fe II or higher ionization 
states are often present. These properties are discussed in, 
e.g. Boroson & Green (1992). In the X-ray band Seyfert 
galaxies are usually well described by a power-law and a 
'soft excess', more emission below ~1 keV than expected 
from an extrapolation of the power-law spectrum observed 
at higher energies. Modelling the soft excess as a black-body 
usually results in a good fit. The soft excess is notable for 
having similar temperatures over a wide range of objects, 
and is the subject of some debate: if it is thermal in origin 
it may be due to a slim disc (e.g Abramowicz et al. 1988), 
Gierlihski & Done (2004) show that soft excess temperature 
in a sample of radio-quiet PG quasars is 0.1 - 0.2 keV over a 
large range of luminosity, and put forward the idea that the 
soft excess is an illusion caused by relativistically blurred 
strong absorption. This paper is based on the idea that it is 
due to photoionized emission blurred by relativistic motion 
in an accretion disc. 

Seyferts are highly variable, in the optical on time-scales of 
weeks to years, and in the X-ray on time-scales of hours or 
even less. NLSls are more variable in the X-ray band than 
typical Seyfert Is, implying they are more compact (Leighly 



1999a). X-rays are generated only in a small energetic area 
of an AGN, i.e. near the central engine. The detection of 
iron lines broadened by relativistic effects in some sources 
shows that the central energy source is matter accreting 
onto a black hole (e.g. Tanaka 1995; Vaughan & Fabian 
2004). 

This concept is extended in relativistically blurred pho- 
toionized disc models, the most recent of which (Ross & 
Fabian 2005) is used in this paper. This model is intended 
to simulate radiation from the inner regions of a black hole 
accretion disc, and it includes atomic physics and general 
relativistic effects. The underlying spectrum is due to a 
slab of optically thick gas of constant density illuminated 
by a power-law, which produces fluorescence lines with 
a continuum of reflected photons. A power-law with the 
spectral index of the illuminating photons is added to 
include light from the illuminating object (e.g. from the 
base of a jet or a hot corona). The power-law component 
dominates the observed spectrum in many NLSls but in PG 
1404+226 the soft excess is the most powerful component 
(Fig. 0. This spectrum is then relativistically blurred as 
calculated for an accretion disc around a Kerr black hole. 
The Laor (1991) line profile (standard in xspec) is used as 
a convolution kernel to do this. 

This model is more physically motivated than the ad hoc 
models usually adopted (i.e. a power-law and a black-body 
to model the soft excess). We show that this new model 
is a better fit to the source than previous models and 
adequately accounts for variability. 



* E-mail: jc@ast.cam.ac.uk 



© 0000 RAS 



2 J. Crummy et al. 




0.5 1 2 5 10 

Energy (keV) 



Figure 1. The ionized disc model as fit to the complete observa- 
tion (see the top of Table l2l for the details of the fit). The dashed 
line is the power-law component, the dotted line is the reflection 
component, and the solid line is their total. 



2 THE SOURCE 

PG 1404+226 is a Narrow Line Seyfert 1 galaxy, with a 
redshift of 0.098. Wang & Lu (2001), using the relation 
with [O III] width, find the central black hole mass to 
be 1.0 x 10 7 M , with an error of ~0.5 dex. During the 
observation we observe a flux of 1.1 x 10 -12 ergs cm -2 s _1 
in the 0.3 - 12.0 keV range, which we convert (using the 
inverse square law with a Hubble parameter of 70 km s" 1 
Mpc -1 ) to a luminosity of 2.3 x 10 43 erg s" 1 . 
The source has previously been observed with the ASCA 
satellite (e.g. Vaughan et al. 1999), and is included in the 
two-part review by Leighly (1999a, b). The ASCA obser- 
vation showed PG 1404+226 to be strongly variable, with 
a clear soft excess. The source showed possible apparent 
absorption around 1 keV, which is discussed in detail in e.g. 
Leighly et al. (1997). PG 1404+226 did not show an iron 
line. As this article will show, the 1 keV feature is a natural 
consequence of an ionized disc. 



3 OBSERVATIONS AND DATA REDUCTION 

PG 1404+226 was observed on-axis with XMM-Newton 
during revolution 0279 (2001-06-17 - 2001-06-18). All EPIC 
cameras were operated in PrimeFullWindow mode with 
the filter Thinl. The RGS instruments were used in the 
HighEventRateWithSES mode. The 0M did not take data. 
During the observation the EPIC MOS took 20.8 ks of data, 
the EPIC pn 18.3 ks and the RGS instruments 21.4 ks. 
The Observation Data Files (ODFs) were reduced in the 
standard way using SAS 5.4.1 to produce event lists. Light 
curves and spectra were created for the sources by using 
a circular, source-centred extraction region of 20 arcsec 
radius for the pn and 30 arcsec for the MOS. The smaller 
than standard extraction radii were chosen to exclude as 
much as possible of the strong background. The background 
light curves and spectra were extracted using regions 
away from any sources and chip-gaps. The presence of 
strong background flaring in the observation was noted and 
compensated for using a Good Time Interval (GTI) file, 



leaving 16.9 ks for the pn and 20.1 ks for the MOS. 
The pn images were also examined for Out-Of-Time (OOT) 
events and pile-up. Any OOT events in PG 1404+226 are 
unobservable due to a low signal to noise ratio, so no OOT 
correction was used. Use of the SAS task epatplot showed 
no pile- up. 

The MOS cameras are not accurately calibrated against the 
pn (Kirsch et al. 2004) , and the pn has a larger effective area 
and spectral range. Therefore all fits were done on pn data 
alone, and checked against the MOS for consistency. The 
MOS gives a flatter spectral index than the pn, but no other 
serious differences were noted. The RGS data were reduced 
and found to be consistent with the other instruments, 
however, the statistics are dominated by the pn and MOS so 
RGS data were not included in any of the fits. The limits in 
accurate calibration were taken as 0.3 keV and 12.0 keV 
for the pn cameras, and the MOS was used over the 0.3 keV 
- 10.0 keV range. The data reduction yielded 7807 (pn) 
source photons (including an estimated 232 background 
photons). 

The spectral analysis was performed using spectra grouped 
so that each bin contained at least 20 source counts, to 
ensure that \ 2 statistics would give reasonable results. 
Response files generated by SAS were used, and the spectra 
were examined using xspec vll.3 (Arnaud 1996). A 
Galactic hydrogen absorption column of 2.14 x 10 20 cm" 2 
was taken from the nh ftool (Dickey & Lockman 1990). Fits 
were also performed allowing for free hydrogen columns at 
the source redshift, no evidence for any extra absorption was 
found. All quoted errors are 90% limits on one parameter. 



4 ANALYSIS 

The data were reduced and spectra were created as de- 
scribed in Section [3] The spectrum is shown in Fig. [5] We 
fit the data with the standard black-body and power-law. 
Adding a red-shifted edge fixed at a source-frame energy 
of 0.87 keV (the O VIII K absorption edge) improves the 
fit (A^ 2 of 18 for 1 degree of freedom). Adding a second 
edge at 0.74 keV (O VII) does not improve the fit. Allowing 
the energy of the first edge to vary improves \ 2 by 4 
for 1 degree of freedom. The results are given in Table 
Q To model cold iron reflection from distant gas (e.g. a 
torus) we added a narrow red-shifted Gaussian at 6.4 keV 
(in the source frame). This improves the fit further, to 
a xt °f 0.916 for 190 degrees of freedom (with all other 
parameters consistent with those reported for the previous 
fit); allowing the width or energy of the line to vary does 
not improve the fit. An F-test gives an 89% probability 
that the addition of the line is valid. Including MOS data 
does not increase the significance, an F-test gives an 84% 
probability that the line is valid when all three EPIC 
observations are fit simultaneously. The use of the F-test 
on hypotheses on the boundary of the parameter space, 
such as testing for the existence of an emission line, is of 
questionable statistical value (Protassov et al. 2002). We 
therefore calibrate the F-test using simulated observations 
(from the xspec 'fakeit' command), by generating data 
from a model without an iron line. In this data we measure 
an iron line at the same or higher significance (measured 
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Figure 2. The pn spectrum of PG 1404+226, showing data 
(crosses) and the ionized disc model fit (line). The lower panel 
shows the residuals to the model, with error bars scaled to a size 
of one. The data has been re-binned for display. 



by an F-test) as detected in our real data 12% of the 
time, therefore we conclude that the F-test is reasonably 
calibrated in this case and that the line is 88% significant. 
The statistics are poor at these energies, so the equivalent 
width of the line is barely constrained; xspec reports the 
equivalent width of the line as 0.5to'g keV. Given that 
the line is just on the boundary of 90% significance, we 
choose not to include the line in any of our further fits. 
We therefore consider the free edge energy model to be our 
best simple fit, with a xt of 0.934 for 191 degrees of freedom. 

We also fit PG 1404+226 using the ionized disc model. 
The free parameters of the model are the iron abundance 
(Fe, in units of solar metallicity), ionization parameter (£, 
the number of incident ionizing photons cm _2 s _1 ), and 
the spectral index of the illuminating continuum (T, this 
is also the index of the power-law component we add). 
The spectrum is blurred with a Laor line profile, which is 
standard in xspec, and has as free parameters the emis- 
sivity index of the disc (Index), the inclination of the disc 
to the line of sight (8, in degrees), the inner radius of the 
disc (Rin, in units of gravitational radii) and outer radius 
(which we fix at 100 gravitational radii as it is not strongly 
constrained). This is intended to reproduce the spectrum 
from a photoionized disc around a maximally-rotating 
(Kerr) black hole. Fits made using relativistic blurring from 
a non-rotating (Schwarzschild) black hole were not found 
to be satisfactory; there were strong residuals below 1 keV 
where the blurring is insufficient to smooth out the emission 
and absorption features. A more rapidly rotating black hole 
allows for stronger blurring. 

The parameters obtained in fitting PG 1404+226 with the 
ionized disc model are given in Table |5J and the model is 
plotted showing the two components in Fig. Q We find no 
evidence for absorption edges with this model. Our best 
ionized disc fit has a xt of 0.892 for 189 degrees of freedom. 



4.1 Variability 

PG 1404+226 varied by a factor of ~3 over the course of the 
observation, as shown in Fig. [3] The Figure shows soft (0.3 - 



1.0 keV) and hard (1.0 - 8.0 keV) count rates. Counts above 
8 keV are not included as the source is extremely weak 
at those energies. The hard count rate does not strongly 
vary whereas the soft count rate increases during the obser- 
vation; this indicates that the spectral shape changes. To 
investigate this change we divide the observation into four 
different periods suggested by the shape of the light-curve, 
as indicated by the vertical lines on the Figure. We fit these 
using our best simple model, testing to see if adding the 
absorption edge (at a fixed energy of 0.9 keV in the source 
frame) improves the fit. The results are given in Table |3] 
Adding an edge results in an improvement to the fit in all 
time periods but the first one. If an edge is included in the 
first period fit XSPEC reports an optical depth of 0.051q os> 
with a negligible change in \ 2 f° r the loss of a degree of 
freedom. To further investigate this we perform simulated 
observations to evaluate whether we would detect an ab- 
sorption edge with a r of 0.9. We measure an optical depth 
> 0.05 in ~99.9% of our simulated observations, therefore 
we are confident that the absence of an absorption edge 
in the first time period is not due to insufficient data quality. 

We also fit the same time periods using the ionized 
disc model, with the results in Table The ionized disc 
model and simple model fit the data about equally well 
when analysing the shorter time periods individually. 
It has been suggested that the variability in some sources 
is due to changes in the normalization of the underlying 
continuum, with a fairly constant reflection component 
from an ionized disc (Fabian et al. 2002, Vaughan & Fabian 
2004). To investigate this we simultaneously fit all four 
time periods, allowing each to have independent reflection 
and power-law component normalizations and tying all the 
other parameters together. We investigate possible changes 
in spectral index of the underlying continuum by repeating 
the previous fit with a free spectral index (F). The results 
from both of these fits are reported in Table [5] We further 
attempt to fit the data with no power-law component at 
all. This did not improve the fits so we do not report the 
results. 

We similarly fit the four periods simultaneously with the 
simple model, tying the black-body temperature, power-law 
index and edge energy together but allowing the normaliza- 
tions and optical depths to vary independently. The results 
are given in Table |3] Our best fit to the data is with the 
ionized disc model, simultaneously fitting all parameters 
except the normalizations and F; this has a xt of 0.918 for 
335 degrees of freedom. 



5 DISCUSSION 
5.1 Simple model 

The simple model is an acceptable fit to the data. This 
model includes absorption at ~0.9 keV, which is likely to 
be due to the O VIII K absorption edge at 0.87 keV. The 
> 0.87 keV energy might imply the absorption originates 
in a slightly relativistic outflow, although the data do not 
constrain this (see Leighly et al. 1997). 

We learn more when we consider spectral variability over 
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Table 1. PG 1404+226 simple model fits. kT is the black-body temperature in keV, Y is the spectral index of the power law component, 
and r is the optical depth of the absorption edge. 



Model 






kT (keV) 


r 


Edge energy (keV) 


r 


xl (d.o.f.) 


power-law - 
power-law - 
power-law - 


- black-body 

- black-body - 

- black-body - 


- edge 
h edge 


1109+ - 0018 
u.±iuy_ Q 0018 

u.±ioo_ 0021 

H44+ - 0019 
u.l±t - .0 019 


i r-o + 0.20 

1 67+ 034 
J -°'-0.31 
i O1+0.19 
1-S1-0.14 


0.87 a 
Q4+0-03 
U - M -0.23 


42+0' 20 
u -* z -0.16 

52+ ' 14 
u - oz -0.13 


1.04(193) 
0.950(192) 
0.934(191) 



a Edge energy fixed. 



Table 2. PG 1404+226 ionized disc model fits, including spectral variation with time. See Section |1] for an explanation of parameters. 



Time (ks) 


Index 




(9 


Fc 


r 


5 


Xl (d.o.f.) 


- 18 a 


Mill 


1 oor+0.362 

l.ZOO_Q Q 


58^24 


5 4+ 2 ' 2 


9 „ 1 7+0.098 
z.oi i _ 123 


87n +430 


0.892(189) 


- 7.5 b 
7.5 - 12 
12 - 14 
14 - 18 


10.0+°° 

2 98+ - 49 
7 q 9 +0.89 
' - a -1.62 
q 07+ - 63 


i ony+0.124 
-0.002 
1 noc + 0.772 

J-.ZOO_Q Q 

1 oqc;+0.332 
i.zoo_ „ 

i 9 o7+0.157 

-0.002 


67 - 3+ 38 

63-2lt 8 4 3 


n 9+12 
Z - -0.60 
1 s+ 8 - 2 
1 - 8 -1.6 

6 - 0+ 27 

2 ll+ i:s 
z - 1J --0.87 


2 40+ - 07 
Z ' 4U -0.25 

9 ona+0192 
2.8U»_ lg0 

9 41q +0.053 
z.<±±y_ 20g 

3 00+° ° 


qcc;+38 

16.8t 47 3 
302.4± 7 ; 2 


0.991(95) 
0.889(92) 
1.12(68) 
0.983(96) 


- 7.5 C 


9 gi+0.09 
a - y -1.64 


, notr+0.065 

i.zoo_ Q 


61-8t? 6 6 


e.o±i:? 


2 401+ - 102 
z.*u±_ 032 


61 n+ 580 

Dlu -220 


0.967(338) 



7.5 - 12 
12 - 14 
14 - 18 



- 7.5 d 9.74±° 3 ; 
7.5 - 12 
12 - 14 
14 - 18 

a The complete observation. b All time periods fit independently. c All time periods fit simultaneously, with all parameters except the 
normalizations of the reflection and power-law components fixed. d As c , but with free spectral index. 



absorption implies the absorbing gas has a high ionization 
parameter; this means it is close to the central black hole 
and the rapid change implied by the variability of the 
feature is plausible. This variation in absorption has been 
previously investigated by Dasgupta et al. (2004), who 
understand it within the model of absorption due to a disc 
wind - as the source brightens more matter is driven off 
in the wind and the absorption intensifies. We have shown 
that the ionized disc model naturally explains this apparent 
increase in absorption. 



5.2 Ionized disc model 

The ionized disc model is the best fit to the entire observa- 
tion. This model has the benefit of being self-consistent and 
not ad hoc. It does not require any additional absorption 
features to give a good fit, the model itself reproduces 
the dip in the spectrum. When the spectral variability of 
the source is investigated using the ionized disc model, it 
accounts for it naturally; the spectral index of the illumi- 
nating powerlaw varies over the course of the observation. 
This fit (at the bottom of Table |5J has a \ 2 20 lower than 
the equivalent simple fit (Table |HJ for two fewer degrees of 
freedom. We conclude that the ionized disc model is the 
most satisfactory explanation of the data. 
The fit parameters give us information about the central 
source. Firstly, it is evident from Table|5|that the inclination 
of the source is quite high; around 60°. This does not fit 



1.2351";? 18 56.6ti 2 2 2 6.1+H 2.163l ;i?J 710± 27 °, 0.918(335) 



9 nno+0.106 

z.ooo_ Q Q51 

9 009+O.IO2 
Z.33Z_q 1Q7 

2 400+ 095 




10 15 

Time (ks) 



Figure 3. The background subtracted light-curve of PG 
1404+226, taken from the pn camera. The circles show soft counts 
from 0.3 - 1.0 keV, and the squares hard counts from 1.0 - 
8.0 keV. The errors are Poisson, and points without associated 
errors have been interpolated over to avoid periods of background 
flaring. The vertical lines show the boundaries between the time 
periods analysed in Tables |21 and 131 



the observation. The absorption appears to 'switch on' 
some way into the observation, as shown in Table [2] The 
absorption increases at the same time as the luminosity of 
the source increases (at ~7.5 ks, see Fig. and we have 
eliminated the possibility that this is a statistical effect. If 
the absorption is due to O VIII, the lack of strong O VII 
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Table 3. PG 1404+226 investigation of spectral variation over time. kT, V and r as in Tabled 'Time' gives the start and end times (in 
ks, counting from the start time) of the section of the observation being fit, as marked on the light-curve given in Fig. [3] 



Time (ks) 


kT (keV) 


r 


Edge energy (keV) 


T 


xl (d.o.f.) 


- 7.5" 
7.5 - 12 
12 - 14 
14 - 18 


1189+ ' 0048 
u.iioa_ 0050 

n 11 7^+0.0054 
u.±i<o_ 0052 

u.nt>«_ 0054 

H03+ 0038 
u.nuo_ 003g 


i cc+0.65 
1 - oo -0.33 
i g 7 +0.33 
± - D '-0.22 
. 74 +0.89 
1 -'^-0.79 
9 n9 +0.80 
z - uz -0.64 


0.9 6 
0.9 6 
0.9 6 
0.9 6 


0.0 C 
n qo+0.36 
u -°°-0.33 
n 77+O.53 
u -' '-0.44 

0.88±g;» 


9871'QQ') 
0.897(95) 
1.14(71) 
0.980(99) 


- 7.5 d 
7.5 - 12 
12 - 14 
14 - 18 


U.11JD_ 0Q24 


1 8^/L+ ai6S 
L854 -0.148 


906+ 1 001 


0.938^ 

645+ ' 407 
u.oto_ 240 

0.986tE!;iJ 


0.971(337) 



a All parameters independent. b Edge energy fixed. c Adding an edge does not improve the fit, so r was fixed at zero. d Investigating 

variation in r; kT, Y and Edge Energy are fit simultaneously. 



very well with the unified model of AGN (Antonucci & 
Miller 1985, Antonucci 1993), in which Seyfert 1 galaxies 
are those visible at low inclinations, with no molecular 
torus in the line of sight. However, PG 1404+226 is fairly 
luminous, and there have been suggestions that brighter 
quasars tend to show less absorption from a torus (Lawrence 
1991, Brandt & Hasinger 2005). The high inclination also 
implies that the brightness of the source is partly due to 
relativistic Doppler beaming from the rotation of the disc. 
We investigated the effect of inclination angle on flux within 
our model, and find that the brightness is decreased by a 
factor of ~5 if the inclination is changed to 0° (all other 
parameters unchanged), and increased by a factor of ~3 
at 90°. The measured luminosity is therefore a reasonable 
order-of-magnitude estimate to the intrinsic luminosity of 
the source. 

Information on the rotation can also be obtained. Note 
that fits were performed assuming a Kerr black hole, and 
that Schwarzschild fits were unacceptable (Section 0J. The 
measured inner disc radius provides further evidence that 
the black hole is rotating - it is at the last stable orbit 
radius for a maximally rotating black hole; well inside the 
plunging region for a non-rotating black hole. The black 
hole rotation is not a free parameter in this analysis so 
we must be careful in drawing conclusions, however it 
seems that a rapidly rotating black hole is likely. It is also 
a possibility that radiation from the plunging region is 
important (see Krolik & Hawley 2002). 

The power-law component of this source is weak compared 
to the reflection component, i.e. the reflection fraction is 
very high. This may be due to light bending effects on the 
primary continuum. If the continuum source is compact and 
close to the black hole almost all the radiation it emits is 
bent onto the disc rather than escaping to the observer; this 
reduces the observed power-law flux while enhancing the 
disc illumination, and therefore the reflection component 
(Miniutti et al. 2003; Miniutti & Fabian 2004). To produce 
a reflection dominated spectrum in this model the primary 
source in PG 1404+226 must be within 3-4 gravitational 
radii of the black hole. Reflected radiation escapes the 
black hole more easily than the illuminating continuum 
as it is emitted by rapidly moving matter and tends to 
be beamed along the plane of the disc to the observer. 
The high inclination therefore also contributes to the high 



reflection fraction. 

This model raises the question of how energy is transported 
through the system. The original source of energy is the 
matter falling into the gravitational potential of the central 
black hole. This heats the electrons that provide the illumi- 
nating continuum, and the energy reflects off the disc and 
escapes. The question of how the in-falling matter heats the 
source of the illuminating continuum is an unsolved one, 
magnetic reconnection is a possibility, as is emission from 
the base of a weak jet, or shocks in a failed jet (Ghisellini, 
Haardt & Matt 2004). In the case of a jet the spin of the 
black hole may also provide energy. 

Given the small disc inner radius we calculate for the 
source, most of the X-ray emission is coming from within a 
few gravitational radii of the black hole, so strong gravity 
effects are expected. This small size also explains the 
rapid variability. If the illuminating continuum is, e.g. UV 
radiation being Compton up-scattered by hot electrons, 
the change in spectral index could be due to a change in 
temperature of the electrons, and the extra flux could be 
due to an increase in emitting area. 



6 CONCLUSIONS 

We have shown that: 

• Physically plausible models based on relativistically 
blurred reflection from photoionized discs fit observations to 
PG 1404+226 better than the typical ad hoc simple model 
of a black-body and power-law. In particular, a 'soft excess' 
at low X-ray energies is a natural consequence of a blurred 
ionized disc; the 'power-law' is largely made up of broad- 
ened iron line emission and the 'soft excess' of other blurred 
lines, plus the Compton reflection component. The ionized 
disc model also leads naturally to apparent absorption fea- 
tures at low X-ray energies. 

• The X-ray emitting accretion disc in PG 1404+226 has 
an inclination of ~ 60°. 

• The central black hole is most likely strongly rotating, 
alternatively, radiation from the plunging region may be im- 
portant. 

• The source has strong and complex spectral variability, 
which is well described within the ionized disc model by a 
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variation in the spectral index of the illuminating contin- 
uum, as well as variation in the flux of the two components. 

• If the alternative simple model is adopted, the spectral 
variation is reasonably fit with a rapid increase in absorption 
at ~0.9 keV coincident with the increase in luminosity of the 
source, plus some variation in the other model parameters. 

The new relativistically blurred photoionized disc model of 
Ross & Fabian (2005) is clearly a valuable tool in inves- 
tigating NLSls, and in sources where an ionized disc fit is 
appropriate information about the central conditions can be 
obtained. The application of these models to more, brighter 
sources (currently in progress) should lead to advances in 
our knowledge of AGN. 
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